Interaction driven real-space condensation 
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We study real-space condensation in a broad class of stochastic mass transport models. We show 
that the steady state of such models has a pair-factorised form which generalizes the standard 
factorized steady states. The condensation in this class of models is driven by interactions which 
give rise to a spatially extended condensate that differs fundamentally from the previously studied 
examples. We present numerical results as well as a theoretical analysis of the condensation transi- 
tion and show that the criterion for condensation is related to the binding-unbinding transition of 
solid-on-solid interfaces. 
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Real-space condensation has been observed in a vari- 
ety of physical contexts such as cluster aggregation [lj , 
jamming in traffic and granular flow [2j, |3j and granular 
clustering |4j. The characteristic feature of these sys- 
tems is the stochastic transport of some conserved quan- 
tity, to be referred to as mass; the condensation transi- 
tion is manifested when above some critical mass den- 
sity a single condensate captures a finite fraction of the 
mass. The condensate corresponds to a dominant cluster 
or a single large jam in these examples. Perhaps more 
surprising realisations of condensation are wealth con- 
densation in macroeconomies ja|, where the condensate 
corresponds to a single individual or enterprise owning 
a finite fraction of the wealth; condensation in growing 
or rewiring networks where a single hub captures a finite 
fraction of the links |6j] and phase separation dynamics in 
one-dimensional driven systems where condensation cor- 
responds to the emergence of a macroscopic domain of 
one phase [7]. 

Mass transport may be modelled in terms of interact- 
ing many-particle systems governed by stochastic dynam- 
ical rules. Generically these systems lack detailed balance 
and thus have nontrivial nonequilibrium steady states. 
Although our understanding of such steady states is still 
at an early stage, a class of models has been determined 
which exhibit a factorised steady state (FSS) Q which 
can be written as a product of factors, one factor for each 
site of the system. This simple form for the steady state 
has afforded an opportunity to study condensation an- 
alytically and has also been used as an approximation 
to more complicated nonequilibrium steady states. The 
conditions under which condensation can occur have been 
determined, leading to conditions on the stochastic mass 
transport rules for condensation to result |3|. One key 
feature of the condensate arising in these models is that 
it forms at a single site. 

In the physical systems of the kind described above, 
generically the stochastic transport rules depend not only 
on the departure site but also on the surrounding envi- 



ronment. In general such models do not have FSS and 
finding their steady states has remained a challenge. 

The purpose of this letter is twofold. First, we intro- 
duce a broad class of mass transport models where the 
transport rules depend on the environment of the depar- 
ture site. These models do not have an FSS, yet we can 
determine their steady states explicitly. The structure of 
the steady state generalises the FSS to a pair-factorised 
steady state (PFSS). Secondly, we find that the nature of 
the condensate in PFSS is strikingly different from that 
of the FSS: unlike in the FSS, the condensate is spatially 
extended. This is due to the short-range correlations in- 
herent in the PFSS, but absent in the FSS. 

We consider a class of mass transport models on a pe- 
riodic chain with sites labelled by i = 1, . . . , L. At each 
site resides a non- negative integer number, to,, of par- 
ticles each of unit mass. We define particle dynamics 
such that a particle hops from site i to i+1 with a rate 
u{mi-i,mi,mi + i) (provided m, > 0), so the total mass 
^2i rrii = M is conserved. These dynamics drive a cur- 
rent of particles through the system. 

If the hop rate is only a function, u(toj), of mass at the 
departure site rm, the model reduces to the zero-range 
process [lOj which has a FSS. Explicitly, the probability 
of a configuration {mi} occurring in the steady state is 
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where f(m) = ^/Y[T=i u (^) for to > 1 and /(0) = 1. 
Thus there is one factor f(rrii) for each site i of the system 
and the delta function ensures that the total mass is M. 
When the hop rates u(m,_i, to^to^+i) depend on all 
three arguments, we propose the PFSS as a natural gen- 
eralization of the FSS which takes the following form: the 
steady state probability of configuration, {mi}, is 
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Thus there is one factor g(mi,mj+i) for each pair of 
neighbouring sites. The normalisation, Zl.m, which 
plays a role analogous to the canonical partition func- 
tion in equilibrium statistical mechanics, is given by 

L L 

zl,m = ^2 n ff ( TO *' m *+i) 5 C^2 mi ~ m ) • ( 3 ) 

{m,}t=l i=l 

Note that in the case g(mi,mi+i) = /(mj), for example, 
the PFSS Eq. © reduces to the FSS form JTJ. 

We first establish that the steady state J5J holds for 
a broad class of mass transport models. We find that if 
(though not only if) the hop rates out of site i factor- 
ize [n|: 

u(toj_i, m*, m»+i) = a(mj_i,mi)/?(mi,mi+i) , (4) 

then the steady state is of the PFSS form (J2J with 
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FIG. 1: Typical steady state configurations for L = 1000 
and J — Uo — 1 (for which p c = 0.51) in (a) condensed 
phase, p = 3, and (b) fluid phase p — 1/4. The inset in (a) 
shows an L 1 / 2 dependence (indicated by the solid line) of the 
condensate width on system size. 
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for to, n > where g(0, 0) = 1, provided a and /3 satisfy 
the constraint 



a(m— 1, n) P(m,n— 1) 
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Furthermore given any form of the weight g(m., ro) one can 
determine the functions a and /3 through the following 
recursions 
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Thus for every choice of g(m, n) there exists a stochas- 
tic mass transport model which will generate the corre- 
sponding PFSS. 

We now focus on a particular model which has a PFSS 
with g(m, n) given by 



g(m, n) = exp 



1 
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One can check from IJ7I4H that the corresponding hop 
rates are 

{exp [—2 J + Uo5 m _i\ , for m < l,n 
exp [2 J + UoS m .i\ , for m>Z,n(9) 

exp [L/o<5m,i] ■ otherwise 

Physically, the rate is low if the mass at the departure 
site is less than the neighboring masses and is high if the 
mass is larger than the neighboring masses. This tends 
to flatten the density profile and generates the effective 
surface tension J in (JHJ) implying short-range correlations 
between the sites. In addition, isolated particles tend 
to hop relatively quickly leading to a preference in the 



steady state weights for vacant sites. This is reflected by 
the on-site attractive potential —Uo8 mit o in (JSJ. 

The model defined by the hop rates @ is guaranteed to 
have a PFSS with g(m, n) in ©. To investigate whether 
the model allows for a condensation transition as the pa- 
rameters J, Uq and the conserved mass density p = M/L 
are varied, we have run Monte Carlo simulations, accord- 
ing to the following prescription. The system is prepared 
in a random, homogeneous initial condition and evolves 
under random sequential update. During each time step 
St a site I is selected randomly and if a particle is present 
it is transferred to the neighbouring site l+l with proba- 
bility u(to;_!, to;, to,; +1 ) St. L such time steps constitute 
a single Monte Carlo step. 

We find that two phases emerge in the steady state 
depending on p and J (where we set Uo = 1). As il- 
lustrated in Fig^ at low density the system resides in 
a fluid phase, in which particles are distributed homo- 
geneously throughout the system. When the density ex- 
ceeds a critical value p c (J), the system is in a condensed 
phase wherein a condensate containing the excess mass 
(p—p c )L coexists with a critical background fluid of mass 
p c L. In contrast to a usual condensate that occupies a 
single site, as for example in an FSS, the condensate here 
extends over many sites. In fact, the condensate extends 
over typically OiL 1 !" 1 ) sites as shown in Fig^(a). 

To locate the phase boundary in the p-J plane we 
computed the single-site probabilities p{m, L) that a site 
contains exactly mass to in the steady state. In the fluid 
phase p(m, L) decays exponentially for large to whereas 
in the condensed phase an additional bump emerges at 
the large to tail of p(m, L) as illustrated in Fig. The 
phase boundary in Fig. |21 is determined by the value of 
J, for fixed p, at which a bump in p(m, L) first appears 
as one increases J . Our theoretical prediction for the 




form of Zl m with respect to M) is given by 



FIG. 2: Single-site distribution p(m, L) for a system of L = 
1000 sites with p — 3 and J — Uq — 1. The inset shows the 
mass at the maximum of the condensate bump as a function 
of system size L, where the crosses are obtained from Monte 
Carlo simulation and the solid line shows an L 1 ' 2 dependence. 
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FIG. 3: Phase diagram for Uo — 1. The crosses are data 
points obtained from Monte Carlo simulation. The solid line 
gives the theoretical prediction 11411 . For J < J c given by (II 12 P 
condensation does not occur at any density. 



phase boundary, presented below, is in excellent agree- 
ment with the numerical results. 

The condensate bump in p(m, L) also has an interest- 
ing scaling behavior with L. We plot in the inset of Fig. 
[5] the mass at the maximum of the bump as a function 
of L and find that it grows as L 1 / 2 . This implies that a 
typical site inside the condensate has mass of order L 1 ! 2 . 
On the other hand since the condensate carries total ex- 
cess mass (p — p c )L this implies that there are typically 
of order L 1 / 2 sites inside the condensate i.e. the spatial 
extent of the condensate is of order L 1 ' 2 . 

We now analyze the general conditions under which 
the steady state © may admit condensation. The grand 
canonical partition function, Zl(p) (the Laplace trans- 
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where the chemical potential p is determined from the 
condition that 



1 din Z L (p) 
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Clearly p(p) is a decreasing function of p for p > |12j . 
If, as p — > 0, the function p{p) — > oo, then a solution 
p > of (fl II) exists for any p > 0. This implies from 
(|TU|l that the single-site mass distribution decays expo- 
nentially for large to signifying a fluid phase and there is 
no condensation. On the other hand, if, as p — > 0, the 
function p{p) approaches a finite value p c , then a solu- 
tion of Qll|l can only be found for p < p c implying that 
the fluid phase exists only for p < p c . When p exceeds 
the critical particle density p c there is no solution to l|ll|) 
implying the onset of condensation wherein the 'excess' 
mass (p — p c )L is carried by the condensate. 

Thus to determine if there is condensation one needs 
to analyze (|10|l and (|ll|l as p — > 0. But, for p = 0, 
(|10|l is precisely the grand canonical partition function 
of a solid-on-solid (SOS) interface model |13, [l4( where 
the interface height at site i is equivalent to the mass 
mi. Since mi > the interface heights are strictly non- 
negative implying that the interface grows on a substrate 
at rrii = 0. Thus p(0) from l|ll|) corresponds to the aver- 
age interface height (rrii) in this SOS model. If p(0) = oo, 
i.e. there is no condensation transition, the interface is 
unbound since its mean height is divergent. On the other 
hand if there is a condensation transition, in which case 
p(0) = p c is finite, the interface is bound with a finite 
mean height p(0). Therefore the criterion for a condensa- 
tion transition is that the corresponding interface should 
be bound. Moreover the critical density p c is given by 
the mean height of the bound interface. 

We now present an exact calculation of the phase di- 
agram in Fig. |21 The mean height of this SOS inter- 
face model can be easily calculated using a standard 
transfer matrix formalism. The partition function 11U|) 
may be written as Zl(p) = Tr [T L (^t)] where the ele- 
ments of the transfer matrix T are g(m' , to). In the large 
L limit only the eigenvector \<fio) of T with the largest 
eigenvalue Ao contributes. The eigenvalue equation reads 
Sm=o5'( m ' ; TO )( m l ( / , a} = A Q (m'|0 Q ). The eigenvectors 
are either extended states (m\(j> a ) ~ e lpm or a bound 
state (to|0 q ) ~ t m where |t| < 1 [l3|- If the spectrum 
contains a bound state then the bound state corresponds 
to the largest eigenvalue. Substituting (m|</>o) = t m into 
the eigenvalue equation for to' > and ml = separately 
yields t = c~ J /(l— e~ u °). For the bound state to exist 
|t| < 1 which implies 

ln{e Uo 



J>Jc = Uo 



1 



(12) 



Therefore, for J < J c the system will not condense 
at any finite density implying p c = oo, whereas for J > 
J c , the system condenses above a finite density p c . The 
density p c is given by the mean height in the bound state 



Pc 
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Using the bound state eigenfunction one finds 



Pc = c^-^-i 



(13) 



(14) 



This prediction is in excellent agreement with the numer- 
ical data as illustrated in Fig. 

We now discuss the condensation transition in a more 
general PFSS where 



g(m, n) — K(\m — n\) exp 



\ (U(m) 



U{n)) 



(15) 



Here K{\m — n\) represents the interaction between near- 
est neighbor masses and —U(m) is an on-site potential. 
If K(x) decays sufficiently rapidly for large x, as in JSJ, 
for condensation to happen one only requires U(m) to be 
positive and localised near m = 0. In this case the con- 
densation is interaction-driven and the existence of the 
condensation transition corresponds to having a bound 
interface. In such cases quite generically the height and 
width of the condensate are expected to scale as L 1 ' 2 . 
This follows from a Brownian excursion argument: the 
localized on-site potential plays no role at sites occupied 
by the condensate — in the absence of the potential, the 
problem can be related to a random walk problem where 
the random walker takes independent steps with length 
drawn from a distribution K(x) 15] . The shape of the 
condensate is determined by a single large loop defined 
by the excursion of the random walker as shown in Fig.^ 
The probability that the walker returns to the origin for 
the first time after N steps scales as N~ 3 / 2 for sufficiently 
rapidly decaying K(x). So, the average number of steps 
until the first return is ff N~ 3/2 N dN ~ L 1/2 (the up- 
per cut-off, L, is determined by the maximum number of 
possible steps). This predicts that the spatial extent of 
the condensate is 0{L x l 2 *). Also, because it is Brownian, 
the typical height of the excursion, and therefore that of 
the condensate, scales as L 1 / 2 . Note that the area under 
the excursion, equivalent to the mass contained in the 
condensate, is O(L), as it should be. 

This interaction-driven condensation is rather differ- 
ent from the type exhibited in an FSS. There the func- 
tion K(x) — 1 and a localized on-site attractive poten- 
tial is no longer capable of driving condensation. Instead 
one requires a specific unbounded potential of the form 
—U(m) ~ 7 In to for large to 10] . Thus in the FSS con- 
densation is potential-driven. 



To summarize, the steady state |J5J extends the class of 
exactly solvable steady states in nonequilibrium statisti- 
cal mechanics. The condensed phase which emerges in a 
PFSS is fundamentally different from that of an FSS. 
In a PFSS the condensation is interaction-driven and 
moreover the condensate extends spatially over OiL 1 ! 2 ) 
sites. The explicit form of the single-site mass distribu- 
tion p(m, L) in the FSS condensed phase has been de- 
termined recently |S( • It remains a challenge to compute 
p(m, L) for the PFSS condensed phase. It would also be 
of interest to study PFSS in higher dimensions. 

Finally we note that the FSS has provided insight into 
number of issues of nonequilibrium statistical physics. 
Although we have focussed here on the issue of conden- 
sation, the generalisation to a PFSS should allow one 
to address other interesting issues such as the role of 
conservation laws [16|. disorder [l?J, boundary-induced 
phenomena [18j and fluctuation theorems (19| . 
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